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Abstract 

The Hamiltonian for a particle constrained to motion near a toroidal helix with 
loops of arbitrary eccentricity is developed. The resulting three dimensional 
Schrodinger equation is reduced to a one dimensional effective equation inclusive 
of curvature effects. A basis set is employed to find low-lying eigenfunctions of 
the helix. Toroidal moments corresponding to the individual eigenfunctions are 
calculated. The dependence of the toroidal moments on the eccentricity of the 
loops is reported. Unlike the classical case, the moments strongly depend on 
the details of loop eccentricity. 

Keywords: toroidal helix, toroidal moment, curvature potential 
1. Introduction 

The majority of work directed towards modeling the metaparticle constituents 
of metamaterials has been performed using classical physics [lj]. The charac- 
teristic length scales of most currently fabricated metaparticles allow for that 
approach to be appropriate and productive. However, it is nearly certain that 
metaparticles will eventually be fabricated on scales at which quantum mechan- 
ical methods will prove necessary to capture their physics with good fidelity 

This paper focuses upon two interesting properties common to many meta- 
particles: they can be approximated as reduced dimensionality systems and they 
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can possess nontrivial topologies. The advent of quasi one and two dimensional 
curved nanostructures has led to situations wherein formalism developed for 
particles constrained to curved manifolds has become of practical importance. 
Specifically, there exits a prescription that allows for degrees of freedom extra- 
neous to the particle's 'motion' on a curve or surface to be shuttled into effective 
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curvature potentials in the Schrodinger equation 

Recently, it was suggested that quantum methods be employed in an ef- 
fort towards understanding toroidal moments induced by currents supported 
on nanoscale metaparticles and the interactions of those moments with time- 
dependent electromagnetic fields 15 1. Because of the theoretical and practical 
interest in toroidal moments Q, 12, 18, 3, 20, 21, 22|, a toroidal helix (TH) 
of adjustable eccentricity has been chosen here to investigate the role of quan- 
tum effects. Being closed, a TH can support current carrying solutions allowing 
for the existence of a toroidal moment [23[. Furthermore, the TH has the ad- 
vantage of having sufficient symmetry to allow for a clean reduction of the full 
Hamiltonian to a one dimensional effective Hamiltonian. 

The goals of this work are threefold. The first is to derive the Hamiltonian for 
a particle in a coordinate system adapted to include points near the coils of a TH 
of arbitrary eccentricity. The next deals with redu cing the full three dimensional 
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25j to arrive at an effective 



Hamiltonian via a well known procedure 
one-dimensional Schrodinger equation. The reduction of dimensionality impels 
the introduction of a curvature potential well known to workers in the field of 
curved manifold quantum mechanics. A basis set consistent with the periodicity 
and symmetry of the system is introduced thereafter. Achieving the first two 
goals and with the basis functions in hand, the spectrum and wave functions 
of the system (which can be used for applications in the external field and/or 
time-dependent case) are found. Finally, toroidal moments corresponding to 
particular eigenstates are determined and their sensitivity to the eccentricity of 
the loops comprising the TH is investigated. 

The remainder of this paper is organized into four sections. Section 2 intro- 
duces a parameterization for an uj turn TH in terms of an azimuthal coordinate 
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tfi. A three dimensional Hamiltonian appropriate to motion near the TH 
follows by attaching a Frenet system to the helix and assigning two coordinates 
Qn-iQb to describe degrees of freedom away from the coil. Section 3 details the 
reduction of to a one-dimensional by standard methods, although per- 
haps unfamiliar to workers in the metamaterial community. As a consequence 
of the reduction, curvature potentials appear. Their presence has been shown 
to be essential in properly describing one dimensional systems that exist in an 



ambient higher dimensional space 



26|. Section 4 presents the basis set used to 



calculate the spectrum, eigenstates and toroidal moments per a given quantum 
state. Those quantities, along with results showing the dependence of TMs on 
eccentricity are given. Section 5 is dedicated to conclusions and some remarks 
concerning future work. 



2. The TH Schrodinger equation 

To arrive at the time independent Schrodinger equation fl^(r) 1 J' = = 
( — i V 2 +V)^f, the Laplacian must be derived from a suitable parameterization 
of the TH geometry. Consider a TH with u> equally spaced circular coils. Let 
R be the distance from the z-axis to a loop center and a the radius of a loop. 
First define 

W{4>) = R + a cos(o;0) (1) 

with <j) the usual cylindrical coordinate azimuthal angle. The circular TH is 
traced out by the Monge form [27 1 

r(<f>) = W((p)p + a sin(w</>)k. (2) 

Generalizing Eq.(2) to coils of arbitrary eccentricity requires only the modifica- 
tion 

r(<f>) = W(4>)p + b sin(a;0)k (3) 

where a,b may be adjusted to yield the coil shape desired (Fig. 1). To avoid 
cluttering the narrative with blocks of equations, the expressions that follow 
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will apply to the circular case only. The expressions for arbitrary a and b are 
given in the appendix. 

A three dimensional neighborhood in the vicinity of the TH is built by 
assigning two coordinates to points near the curve along unit vectors orthogonal 
to the curve's tangent and to each other. The Frenet-Serret equations 27 1 



provide such an orthonormal coordinate system known as a Frenet trihedron. 
The unit tangent to any point on a curve traced by r(<^>) is 



dr(<P) 



dr(<t>) 



from which the Frenet trihedron can be constructed via the relations 



dr(4>) 



dN 



dr(<j>) 



dB 



k(»N 
(-«(0)f + r(0)B) 

dv{4>) 



(4) 

(5) 

(6) 
(7) 



where the curvature and torsion of the space curve r(0) are indicated by k(<P) and 
t(<P) respectively (where again, detailed forms for the expressions in Eqs. (4-7) 
appear in the appendix). Points near the TH are located via two perpendicular 
displacements qjvN and q^B. The TH position vector may now be written 



x(0, q N , q B ) = r(<f>) + g^N + q B B. 



(8) 



It should be noted that Eq.(8) dehnes a Cartesian region about a curve traced 
by r ((/)). While it is certainly possible to construct a hnitc tubular neighbor- 
hood about r (</>), the coordinate ambiguity of the azimuthal angle as the radial 
distance approaches zero causes the limiting procedure to become complicated. 
Additionally, the separability of the Schrodinger equation into tangential and 
normal variables is lost, and with it any real advantage in using the reduced 
Hamiltonian. 

The covariant metric tensor elements Qij can be read off of the quadratic 



form 



28] 



dx ■ (ix = gijdq l dq J 



(9) 
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where in what follows the ordering convention is (g 1 ,^ 2 ,^ 3 ) = (<f>, q^, q B )- The 
Laplacian is 

with g — det(gij) and g^ the contravariant components of the metric tensor. 
Before presenting explicit forms for g^ and g lj , it is useful to define 

dv{4>) 



d(j> 



[(auf + W{<f>) 2 ] 1 ' 2 



and 



G((f>,q N ) = 1 - q N n(4>) 
after which the covariant metric may be written 



(11) 



(12) 



9ij 



(f(mG^qN) 2 +r^n q % + q %)} -T{<t>)q B f{<t>) r(<f>)q N f(^ 
-T(<f>)q B f(cj>) 1 

r(<p)q N f(<P) 1 ) 



The contravariant form of the metric is obtained straightforwardly; 



(13) 



f(<j>) 2 G(<j>,q N ) 2 



1 r(4>)q B f^) -T^)q N f(4>) \ 

T(4)q B f(4>) m 2 [G(<P,qN) 2 +r(cP) 2 q%} -T(cj>) 2 q N q B f(<j>) 2 

\2„._„„ ft A\2 flAXt\n(A „..12 , Jl( 



\-T(<}>)qNf(ct>) -T((f>) 2 q N q B f(<f>f 



/(0) 2 [G(^ 9Af ) 2 + r 2 (0) 9 2 r ]y 
(14) 



It is easy to show that 



V5 = /W(1-?W^))- 



The Laplacian found by directly evaluating Eq.(lO) is complicated by the exis- 
tence of cross terms arising from d 2 /dq l dq : > , (i ^ j), operations. However, all of 
those terms arc multiplied by the distance parameters qw and q B such that in 
the limit qN,qB — > they vanish independently of the derivative operators that 
follow them. Taking this limit now (it will be taken again later post operation 
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of the Qn,b derivatives) leads to a more convenient starting point for develop- 
ing the reduced Hamiltonian in the ensuing section. Physically, the limiting 
procedure is effected by external mechanical or electrical constraints; mathe- 
matically, they are added by hand into the Schrodinger equation as potentials 
V n (q) normal to the lower dimensionality base manifold. Their detailed forms 
are not important. Previous work has shown that even for finite thicknesses, 
degrees of freedom extraneous to those of the base manifold do not mix with 
the latter in the sense that their wave functions decouple |26j. Here, for the 
sake of dcfmitcness, hard wall potentials are assumed for V n (q lss ) and V n (q B ) in 
this and the next section. With this discussion in mind, H% may be written as 
(with h = m = 1) 



If 1 d 2 f{4>) d d d 2 d 2 

k ) 

(15) 

Note that the at this stage is still not separable. The procedure for rendering 
H% separable and arriving at a simpler effective Hamiltonian is given in the 
following section. 

3. Constructing the effective Hamiltonian 

As the particle is constrained to the toroidal helix, its wave function will 
decouple into tangent and normal functions (the subscripts t and n denote 
tangent and normal respectively) 

*(0,97V,te) Xt{4>)Xn(qN)Xn{qB) (16) 

and G(c/f>, qw) will approach unity. The normalization condition 

/ \^(c/ ) ,q N ,q B )\ 2 G( ( j > ,q N )f^)dcbdq N dq B = l (17) 
Jo 

becomes 



2tt 

IXtWOI 2 \Xn(q N )\ 2 \Xn(q B )\ 2 f(<P)dcPdq N dq B = 1. (18) 
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The norm must be conserved in the decoupled limit 6] , which implies 

\^{4>,q N ,q B )\ 2 G{^q N ) = \xt^)\ 2 \Xn(q N )\ 2 \Xn(q B )\ 2 ■ (19) 

The wave function ^(4>,q N ,q B ) is now related to Xt(4>)Xn{qN)Xn(q B ) by 

^{<p,qN,q B ) ^ Xt{<t>)Xn{qN)Xn{q B )G- 1/2 {(j),q N ). (20) 

Applying to ^((/), g^r ,q B ) and taking the limit as qw, q B ~ > post all deriva- 
tive operations yields the result 

tt3 l ( 1 d 2 /'(</)) a i 2 ,,, a 2 a 2 \ T _, , T , . , 

H - = -AlWW~lWd4> + ^ W+ dqJ + dql) +VniqN) + Vn{qB) - 

(21) 

Distributing the energy between the (cj), q^, q B ) degrees of freedom by allowing 
E = E,f, + En + E B , leads to the decoupled system 

' / '"' ) 3 +±K 2 (4,))xt(t)=E <t>Xt (4>) (22) 



2V/(0) 2 50 2 f^fdcj, 4 
1 d 2 Xn(qN) 



2 dq 2 N 
1 d 2 Xn(qB) 



Vn{qN)Xn(qN) = E N Xn (<?AT ) (23) 



+ K(<7b)x«(9b) = E B x n {q B ). (24) 



2 a?! 

Since ^(gAr) and V a ((J , _b) are the confining potentials effecting the qN,q B — > 
constraint, qjy and q B can be considered spectator variables and only the 
^-dependent part of the Hamiltonian indicated in Eq.(21) is nontrivial. The 
Hamiltonian in one dimension Hr, is written 



with 



V c {4>) = -i K 2 (0) (26) 



the curvature potential. The curvature potential V c (<p) emerges as an artifact of 
embedding the particle's one dimensional path of motion in the ambient three 
dimensional space. The explicit form of the curvature potential in Eq.(26) can 
be determined from 

^) = [P 1 (0) 2 +P 2 (<£) 2 ] 1/2 (27) 
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where 



and 



auS 1 + W(<p) cos (oj(j 



W) = 



sin(w^) 

~1W 



1 



2 -| 



(28) 
(29) 



Explicit forms of the tangent, normal, and binormal vectors, along with other 
vectors and functions for the circular and elliptic helices are given in the ap- 
pendix. 

A plot of V c (<p) for some representative values of a, b with uj = 4 appears in 
Fig. 2. Note that the circular case values are negligible in magnitude compared 
to the eccentric cases, and when a > b, V c {4>) is substantially larger than for the 
converse. For larger ratios of a to b, V c (<p) can be orders of magnitude larger 
than indicated in the figure. 

It is worth stating that instead of parameterizing the TH with 0, it would 
also be possible to employ an arc length scheme where an arc length parameter A 
is determined from A = f Q /(</>') dip' . However, to include the curvature potential 
as a function of A, it would be necessary to find <p(\) along the curve. While 
this could be accomplished numerically, using the azimuthal angle is somewhat 
better suited to incorporating external fields 29 



3- 



4. Computational methods and results 

If the TH is small enough to require a quantum mechanical description, 
the (^-dependent part of its wave function must obey Bloch's theorem (the t- 
subscript will be dropped hereafter) 

Xk(<P). (30) 



A standard choice is 



31] 











= exp 


ik — 






lj 



Xk(<j>) = exp(ik(p)uk(<f>) (31) 

where Uk(<p + 2tt/uj) = Uk(<p) is satisfied. Single valuedness requires the Bloch 
index k = p = integer. A convenient choice for u p (<p) basis elements is 

UpGM) =exp[t j 00]. (32) 
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The requirement indicated in Eq.(30) yields 

(3 = um, n = integer. (33) 
From the above considerations, a suitable basis set for the TH is 

x pa {<f) = explipcp] ]T CTexp(ina;0). (34) 

n 

The Bloch form introduces sub-states (sub-bands in the case of a continuous 
rather than discrete index) for each p value which would not be present if the 
TH were treated as a ring of length L. The C£f are the expansion coefficients 
for cc-th sub-state of a given p value. In this work, it was found that a five-state 
expansion proved sufficient to yield basis size independent results for the lower 
p sub-states. For u turns, values of p consistent with the Bloch theorem, p < u, 
are used. For clarity, only p > are discussed. 

A disadvantage of directly adopting the expression given by Eq. (34) is that 
the basis functions are not orthogonal over the integration measure f(<f>)d<f>. A 
more natural basis set is given by a re-scaled form of Eq. (34) 

exp[ip<p] 



^Ecrexp(^). (35) 



With basis function orthogonality preserved on the right hand side of the Hamil- 
tonian, eigenvalues and eigenvectors are calculated by diagonalizing the matrix 
comprising the elements 



(p + umf + 5f'(d>r f"(cP) 



lo ^ fW 4 /(</,)<> 2/(0)3 (3g) 

-2^ + ^)^] dep. 

Once the eigenstates are found, the current in general is calculated with (now 
with units) 

j((f>,q B ,q N ) - — Im[^*{(j),q B ,q N )V^{(j),q B ,q N )\. (37) 
m e 

The current density given by Eq.(37) is inclusive of cross-sectional degrees of 
freedom and yields a current passing through a rectangular area with unit nor- 
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mal T. However, in keeping with the intent of this work, the limit of infinites- 
imal thickness is assumed (or equivalently, the qb , Qn degrees of freedom are 
integrated out) leading to the current expression for the pa-th state 



j pQ (0,O,O) = ^/m 



(x pQ (</>))* dx^tt) 



(38) 



where the form of the reduced gradient operator is obvious. The quantum 
mechanical current that stems from Eqs.(35) and (38) becomes 

n<t>) „ s . 



(p + ujn) 



2/(0) 



3 sin[a;(ri— m)(f, 
(39) 



When V c ((f>) is included in the Hamiltonian the C^f are modified, causing the 
current to become inclusive of curvature effects. This current is then used to 
calculate the toroidal moments according to [32| 

f27T 



± M 



1 

10 J 



[(j^tf, 0, 0) • r)r - 2r 2 f a (0, 0, 0) ] /(<£)#. 



(40) 



Equation (40) allows calculation of quantum mechanical toroidal moments of 
ground and excited states for each Block index p. For a macroscopic thin wire 
where jdr — ¥ Idr is applicable, the toroidal moment for each p reduces to the 
classical result 



rpp _ J_ 
10 



,dr , „ 7 dr 

^■ V)r - 2r d4> 



f{4>)d<t>. 



For circular TH, Eq.(41) yields 



= iiu I a 2 R c 



(41) 



(42) 



and for the elliptic TH, 



rpp 

M 



■KuIabR 



(43) 



As a means of comparison, the current for the p state without curvature effects 
(i.e. a free particle on a given uj turn helix) is easily determined to be 

2nq e hp 



I = 



„L 2 



(44) 
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where the total length of the TH, L, is calculated using 

Jo 

The formalism described in this section was employed to calculate the eigen- 
values and eigenstates expressed in terms of the C£f for several to and p values. 
To get a sense of the modifications arising from V c (<fi), the eigenvalues and am- 
plitudes for a six-turn eccentric helix in a p = 1 state are listed without (Table 
1) and with (Table 2) the curvature potential being present. The eigenvalue 
shifts reflect that V c {4>) is always attractive as shown in Fig. (2), and capable of 
causing amplitude shifts. The reader will note there is no table indicating the 
shifts for the circular case; the effects are negligible and essentially independent 
of the coil radius a. 

With the CP" amplitudes in hand, Eq.(39) can be used to find the j pa ((f>, 0, 0) 
necessary for computing TMs. To set a baseline for understanding the effect of 
including V c ((j)), the curvature potential was shut off and j pa ((f), 0, 0) determined 
for many cases. In Fig. (3), representative results are given for w = 4, p = 1. 
As anticipated, the lowest energy states yield very steady currents; oscillations 
begin to manifest in the higher sub-band energy states. Turning the potential 
on produces very little change in the currents; in Fig. (4), it becomes clear that 
V c ((j)) does little, which is consistent with its small amplitude indicated in Fig. 
(2). 

When eccentricity is introduced by setting a = .75 and b = .25, the results 
are less trivial. The results displayed in Figs. (5) and (6) are representative 
of a general trend observed throughout values of cj,p. The curvature potential 
suppresses the current in every sub-band by a discernable fraction. Similar 
behavior is observed when a = .25 and b = .75 as shown in Fig. (7) and (8), but 
note that the magnitudes are substantially different from the converse values 
of a and b. The Bloch form of the wave function, independent of the presence 
of V c {(j)) (which again lessens the magnitude of the currents), the Bloch form 
of the wave function and the ab dependence of the Laplacian are sufficient to 
cause asymmetries in the currents. 
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Toroidal moment results for lj = 4 are shown in Table 3 for several p-states 
and their corresponding sub-states. For a given p value, the lowest energy state 
in Table 3 agrees very well with the value obtained if the current given by 
Eq.(44) were used, although V c (<j>) can shift the ordering of states in a way to 
reorder the ground state moment as seen for the p = 3 states. In isolation, 
this is relatively unimportant. However, in a broader context where the natural 
temperature scale of a lOOOA helix is a few \xK, thermodynamic averages of the 
type 

(T M ) = J2 T M (E n )exp[-E n /T] (45) 

will necessitate accounting for proper ordering. 

The modifications to the TMs for the upright (b > a) coil situation are 
generally minimal with exceptions only for p — 2. The flattened coil (a > b) 
results in Table 3 show a much stronger variation in TM values, consistent with 
the much larger strength of V c (4>) for a > b relative to the converse. A sense of 
the dependence of TMs on co can be gleaned from Table 4 where now u = 8. 
Increased variation is seen for both eccentricities, but the flattened coil case 
demonstrates appreciable deviation from the classical expression. 

5. Conclusions 

In this work a prescription to include curvature at the nanoscale for particles 
constrained to toroidal helices was presented, which the authors applied toward 
a quantum mechanical calculation of toroidal moments. It is worth emphasizing 
that the curvature inclusive reduced dimensionality Schrodinger equation devel- 
oped here is driven by an interest in having more tractable, effective models for 
nanomaterials, and is done with the aim of eventually confronting experimental 
data rather than as a purely theoretical exercise. In that context, the choice to 
consider helices was driven by their capability of producing toroidal moments, 
which are currently of both theoretical and practical interest. The curvature 
potential for the helix was derived and shown to be the dominant part of the 
Hamiltonian for lower energy eigenstates of eccentric helices. An intriguing 
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result that arose here was a demonstrated ab asymmetry in several states of 
the quantum mechanically calculated Tj|", an asymmetry not exhibited in the 
classical expression of Eq.(43). 

The array of results given in this work was limited to relatively small values of 
u> and to less severe eccentricity because of numerical limitations on evaluating 
integrals of the type shown in Eq.(36). The extension to larger values of u 
were considered (at least currently) outside the scope of what the authors were 
attempting to accomplish. However, preliminary work gives some indication 
that Mathematica is capable of performing the necessary integrals, albeit with 
increased time expense. It would be of interest to investigate more extreme cases 
of eccentricity and loop number given the enhancement of moments already 
evidenced by larger cj. 

Tailoring the response of toroidal helices to electromagnetic radiation by 
fabricating objects with curvature as a free parameter is still well outside the 
reach of current fabrication methods. However, the formalism and basis set 
established here may serve as means for further investigation of the Tm • d~E/dt 
interactions relevant to the coupling of toroidal moments to electromagnetic 
fields. The extension of the methods here to cases where external vector po- 
tentials are present may be naturally developed from work already done for 
tori immersed in arbitrary magnetic fields [29| and is ongoing with an aim to 
understanding persistent current effects. 

Finally, debate as to whether curvature effects are relevant to, and how they 
arc manifested in, topologically novel nanostructures may eventually be settled 
by examining systems akin to toroidal helices. By opting to either include 
or exclude curvature potentials in modeling routines, it may prove true that 
sensitive quantities like toroidal moments will provide a clear signature as to the 
influence of V c . Work such as that done in this paper may hopefully contribute 
to a resolution to the question of how to properly incorporate twists and turns 
in the quantum mechanical description of bent nanostructures. 
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Appendix 

This appendix presents a more complete set of formulae for the circular TH 
as well as a corresponding set for the elliptical case. 

A.l The circular case 



9 = —sm(uj^))p + cos(cj0)k 
n = cos(u)(f>)p + sin(w0)k 



e 2 



W((f>)0 - aucj) 
W) 



auj 2 + W((j>)cos(u!(l>) 



N : 

B = 



f{4>) 2 
sm(u><j)) ( acu \ 

~JWl + \1W)j . 

/c(0 = (i?M+P 2 W /2 
auj9 + W{4>)4> 

7¥) 

(P2(^)e 2 + Pi(^)n) 



T = 
1 

1 



(_P l( 0)e 2 +P 2 (0) n ) 



(A-l 
(A-2 
(A-3 
(A-4 

(A-5 

(A-6 
(A-7 
(A-8 

(A-9 

(A-10 



A. 2 The elliptic case 

With W(4>) = R + a cos(w^), the equation of the elliptic toroidal helix is 

r(<f>) = W{4>)p + b sin(w^>)k. 

The extension of the results of Sec. A.l to the elliptic case is straightforward: 

Y' 2 (A-ll) 



P((j)) = [ a 2 sin 2 (u;</>) + b 2 cos^(w</ 
1 



9 R — 



[—a sm(w<f>)p + b cos(w0)k] 



(A-12) 
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n£ = [b cos(w0)p + a sin(u;</>)k] (A-13) 

62 = 7w (A " 14) 

ftf) = {P{4>) 2 lo 2 + WW 2 ) 1 ' 2 (A-15) 
6 a^ 2 + VK(</>)cos(^) 

PlW -"pw iw (A " 16) 

sin(w0)/_a_ u 2 iy(0)(a 2 - 6 2 )cos(a;^) + P(c/)) 2 auj 2 \ 

«(0) = (P 1 (0) 2 +P 2 (^) 2 ) 1 / 2 (A-18) 

- _ P(4>)u e E + w(4>)j> 



(A-19) 



N = —y-z (P2 (0)e 2 + Pi (<f>)h E ) (A-20) 
«(<?) 

B = ^y(-Pi(^)e 2 + P 2 (0)n £ ) (A-21) 
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Fig. 1: A toroidal helix where R is the distance from the center of the TH to a 
center of a loop of the TH. The parameter a is the TH's maximum perpendicular 
horizontal distance from a concentric cylinder of radius R. The parameter b is 
the TH's maximum vertical distance from the x-y plane. 
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Fig. 2: The curvature potential V c (4>) in units of h 2 / (m e R 2 ) for the case of the 
circular TH with R = 1, ui = 4, a = b = 0.5 and two elliptic TH cases: R = 1, 
uj = 4, a = 0.25, b = 0.75 and R = 1, u> = 4, a = 0.25, 6 = 0.75. 
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Fig. 3: j P "(0,O, 0) = in units of q e U/(m e R 2 ) for five eigenstates of the 
circular TH configuration lj = 4, a = 0.5, b = 0.5, R = 1, p = 1 where V c ((f>) is 
neglected. 
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Fig. 4: j pa ((f> 7 0, 0) = I((f>) in units of q e h/(m e R 2 ) for five eigenstates of the 
circular TH configuration co = 4, a = 0.5, b = 0.5, R = 1, p = 1 with Vc(</>)- 
Curvature effects on the current are small in the circular case. 
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Fig. 5: j p "(</>, 0, 0) = I{<p) corresponding to the ground state, and second and 
fourth excited states in units of q e h/{m e R 2 ), for the elliptic TH configuration 
w = 4, a = 0.75, b = 0.25, R = 1, p = 1, without the curvature potential V c (</>). 
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Fig. 6: j pa {4> 1 0, 0) = I{<p) corresponding to the ground state, and second and 
fourth excited states in units of q e h/(m e R 2 ), for the elliptic TH configuration 
w = 4, a = 0.75, b = 0.25, R = 1, p = 1 inclusive of the curvature potential 
T4(0). Inclusion of the curvature potential causes a reduction in amplitude for 
each j. 
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Fig. 7: j pa (0,0,0) = I{<p) corresponding to the ground state, and second and 
fourth excited states in units of q e h/(m e R 2 ) for the elliptic TH configuration 
lu = 4, a = 0.25, b = 0.75, R = 1, p — 1 without curvature potential. 
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Fig. 8: j p "(i/),0, 0) = I{4>) corresponding to the ground state, and second 
and fourth excited states in units of q e h/(m e R 2 ) for the configuration w = 4, 
a = 0.25, b = 0.75, i? = 1, p = 1 inclusive of the curvature potential. When the 
curvature potential is included in the calculation of /(</>), the current is reduced 
in amplitude. 
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Eo 


El 


E2 


E3 


Ei 




0.0724 


1.6369 


3.0045 


7.7907 


10.9186 


m 


£,(0) 


ci 1} 


ci 2) 


ci 3) 




-2 


-0.1055 


-0.1648 


-0.0607 


-0.9761 


0.0722 


-1 


0.0585 


-0.9762 


-0.1236 


0.1631 


-0.0428 





0.9822 


0.0315 


0.0291 


-0.1020 


0.1520 


1 


0.0556 


0.1374 


-0.9702 


0.0171 


-0.1907 


2 


-0.1331 


-0.0087 


-0.1970 


0.0996 


0.9662 



Table 1: Eigenvalues and amplitudes Cm for an us=6, a=0.75, 6=0.25, R=l, 
p=l elliptic TH neglecting curvature effects. 





Eq 


Ei 


E2 


^3 


Ei 




-1.4739 


-0.0442 


2.1393 


5.7258 


9.2713 


m 






r (2) 






-2 


0.0606 


-0.0310 


0.1231 


-0.9410 


0.3077 


-1 


-0.3794 


-0.7975 


0.4620 


0.0389 


-0.0713 





0.8927 


-0.4441 


-0.0424 


0.0579 


-0.0266 


1 


-0.2353 


-0.4071 


-0.8712 


-0.0772 


0.1179 


2 


-0.0062 


0.0118 


-0.1027 


-0.3220 


-0.9411 



Table 2: Eigenvalues and amplitudes Cm for an cu—6, a=0.75, 6=0.25, R=l, 
p=l elliptic TH including curvature effects. 



26 



Table 3: Toroidal moments for two configurations with oj=4. TM's with and 
without curvature effects, and classical calculation for each case. 



w 


a 


b 






4 


0.25 


0.75 






P 


TM 


TM wlV c (4>) 


ratio 


Classical TM 


1 


-0.0334 


-0.0317 


1.0545 


-0.0332 


1 


0.0895 


0.0718 


1.2457 




1 


-0.1478 


-0.1308 


1.1293 




1 


0.1901 


0.1837 


1.0347 




1 


-0.2401 


-0.2347 


1.0230 




2 


-0.0669 


-0.0551 


1.2129 


-0.0664 


2 


0.0600 


0.0524 


1.1452 




2 


-0.0088 


0.0193 


-0.4562 






-0.0031 


-0.0344 


0.0896 




2 


-0.2646 


-0.2655 


0.9967 




3 


0.0288 


0.0305 


0.9430 


-0.0995 


3 


-0.0993 


-0.0798 


1.2441 




3 


0.1380 


0.1203 


1.1479 




3 


-0.2038 


-0.2052 


0.9930 




3 


-0.2888 


-0.2908 


0.9931 





a 


a 


b 






4 


0.75 


0.25 






P 


TM 


TM w/V c (<j>) 


ratio 


Classical TM 


1 


-0.0317 


-0.0068 


4.6359 


-0.0319 


1 


0.1625 


0.0627 


2.5936 




1 


-0.2697 


-0.1743 


1.5468 




1 


0.1505 


0.1470 


1.0240 




1 


-0.1694 


-0.1863 


0.9096 




2 


-0.0561 


0.0000 




-0.0638 


2 


0.1067 


0.0157 


6.7927 




2 


-0.3206 


-0.1161 


2.7623 




2 


0.1298 


-0.0081 


-15.9737 




2 


-0.1755 


-0.2072 


0.8472 




3 


0.0664 


0.0063 


10.5665 


-0.0957 


3 


-0.0933 


-0.0346 


2.6957 




3 


0.1211 


0.1092 


1.1090 




3 


-0.3893 


-0.3413 


1.1407 




3 


-0.1784 


-0.2131 


0.8372 
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Table 4: Toroidal moments for two configurations with w=8. TM's with and 
without curvature effects, and classical calculation for each case. 





a 


b 






8 


0.25 


0.75 






P 


TM 


TM w/VcM 


ratio 


Classical TM 


1 


-0.0190 


-0.0166 


1.1485 


-0.0195 


1 


0.1218 


0.0441 


2.7619 




1 


-0.1600 


-0.0822 


1.9477 




1 


0.2787 


0.1938 


1.4383 




1 


-0.3175 


-0.2351 


1.3501 




2 


-0.0386 


-0.0332 


1.1615 


-0.0390 


2 


0.1123 


0.0664 


1.6917 




2 


-0.1890 


-0.1430 


1.3217 




2 


0.2635 


0.2293 


1.1490 




2 


-0.3401 


-0.3113 


1.0923 




3 


-0.0582 


-0.0490 


1.1888 


-0.0584 


3 


0.0952 


0.0699 


1.3622 




3 


-0.2104 


-0.1853 


1.1354 




3 


0.2454 


0.2238 


1.0966 




3 


-0.3598 


-0.3472 


1.0363 







a 


6 






8 


0.75 


0.25 






P 


TM 


TM w/140) 


ratio 


Classical TM 


1 


-0.0189 


-0.0096 


1.9789 


-0.0192 


1 


0.2454 


0.0560 


4.3789 




1 


-0.3224 


-0.1299 


2.4826 




1 


0.3520 


0.2586 


1.3612 




1 


-0.3971 


-0.3162 


1.2558 




2 


-0.0378 


-0.0157 


2.4149 


-0.0383 


2 


0.2267 


0.0757 


2.9954 




2 


-0.3806 


-0.2231 


1.7056 




2 


0.3382 


0.2988 


1.1315 




2 


-0.4284 


-0.4176 


1.0257 




3 


-0.0565 


-0.0138 


4.1053 


-0.0575 


3 


0.1951 


0.0614 


3.1800 




3 


-0.4221 


-0.2692 


1.5680 




3 


0.3116 


0.2627 


1.1858 




3 


-0.4509 


-0.4640 


0.9718 
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